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We present the results of ~2.5 years AGILE observations of PSRB1509-58and of the same interval 
of Fermi-LAT observations. The modulation significance of AGILE light-curve above 30 MeV is at 
a 5(T confidence level and the light-curve shows a broad asymmetric first peak reaching its maximum 
0.39 ± 0.02 cycles after the radio peak plus a second peak at 0.94 ± 0.03. The gamma-ray spectral 
energy distribution of pulsed flux is well described by a power-law (photon index a = f .87 ± 0.09) 
with a remarkable cutoff at iJc = 8f ± 20 MeV, representing the softest spectrum observed among 
7-ray pulsars so far. The unusual soft break in the spectrum of PSRBf509-58 has been interpreted 
in the framework of polar cap models as a signature of the exotic photon splitting process in the 
strong magnetic field of this pulsar. In the case of an outer-gap scenario, or the two pole caustic 
model, better constraints on the geometry of the emission would be needed from the radio band in 
order to establish whether the conditions required by the models to reproduce AGILE light-curves 
and spectra match the polarization measurements. 



I. INTRODUCTION 

PSRB1509-58 was discovered as an X-ray pulsar 
with the Einstein satellite and soon also detected at 
radio frequencies [26j . with a derived distance sup- 
porting the association with the SNR MSH 15-52 
{d ^ 5.2 kpc). With a period P ~ 150 ms and a pe- 
riod derivative P ~ 1.53 x 10~^^s s""'^, assuming the 
standard dipole vacuum model, the estimated spin- 
down age for this pulsar is 1570 years and its inferred 
surface magnetic field is one of the highest observed 
for an ordinary radio pulsar: B — 3.1 x 10^'^ G, as 
calculated at the pole. Its rotational energy loss rate 
is £; = 1.8 X 1037 erg/s. 

The young age and the high rotational energy 
loss rate made this pulsar a promising target for 
the gamma-ray satellites. In fact, the instruments 
on board of the Gompton Gamma-Ray Observatory 
(GGRO) observed its pulsation at low gamma-ray 
energies, but it was not detected with high signifi- 
cance by the Energetic Gamma-Ray Experiment Tele- 
scope {EGRET), the instrument operating at the en- 
ergies from 30 MeV to 30 GeV. This was remarkable, 
since all other known gamma-ray pulsars show spec- 
tral turnovers well above 100 MeV (e.g. [11]). Hard- 
ing et al. ([131) suggested that the break in the spec- 
trum could be interpreted as due to inhibition of the 
pair-production caused by the photon-splitting phe- 
nomenon [3]. The photon splitting appears, in the 
frame of the polar cap models, in relation with a very 
high magnetic field. An alternative explanation is pro- 
posed by ^501 using a three dimensional outer gap 
model. They propose that the gamma-ray emission 
is produced by synchrotron-self Gompton radiation 



above the outer gap. 

The Itahan satellite AGILE [41*1 obtained the first 
detection of PSRB1509-58 in the EGREThand (Pel- 
lizzoni et al. 2009b) confirming the occurrence of a 
spectral break. Here we summarize the results of 
a ^ 2.5 yr monitoring campaign of PSRB1509-58 
with AGILE, improving counts statistics, and there- 
fore light curve characterization, with respect to ear- 
lier AGILE observations. More details on this analy- 
sis can be found in [sBJ. With these observations the 
spectral energy distribution (SED) at i? < 300 MeV, 
where the remarkable spectral turnover is observed, 
can be assessed. 



II. OBSERVATIONS, DATA ANALYSES AND 
RESULTS 



AGILE devoted a large amount of observing time 
to the region of PSRB1509-58. For details on AGILE 
observing strategy, timing calibration and gamma-ray 
pulsars analysis the reader can refer to |3^. A 
total exposure of 3.8 x 10^ cm^ s {E > 100 MeV) was 
obtained during the 2.5 yr period of observations (July 
2007 - October 2009) which, combined with AGILE 
effective area, gives our observations a good photon 
harvest from this pulsar. 

Simultaneous radio observations of PSRB1509-58 
with the Parkes radiotelescope in Australia are ongo- 
ing since the epoch of AGILE^s launch. Strong timing 
noise was present and it was accounted for using the 
fitwaves technique developed in the framework of the 
TEMP02 radio timing software [H, Using the 
radio ephemcris provided by the Parkes telescope, we 
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FIG. 1: Phase-aligned gamma-ray light curves of 
PSRB1509-58 with radio peak at phase 0. From the top: 
AGILE > 100 MeV, 20 bins, 7.5 ms resolution; AGILE 
< 100 MeV, 10 bins, 15 ms resolution; COMPTEL 10- 
30 MeV and COMPTEL 0.75-30 MeV (from [H). 



performed the folding of the gamma-ray light curve 
including the wave terms [33| . An optimized analy- 
sis followed, aimed at cross-checking and maximiza- 
tion of the significance of the detection, including an 
energy- dependent events extraction angle around the 
source position based on the instrument point-spread- 
function (PSF). The chi-squared (x^)-test applied to 
the 10 bin light curve at i? > 30 MeV gave a detec- 
tion significance of ct = 4.8. The unbinned Z^-test 
gave a significance of ct = 5.0 with n = 2 harmon- 
ics. The difference between the radio and gamma- 



ray ephemerides was APradio,' 



10 ^ s, at a level 



lower than the error in the parameter, showing perfect 
agreement among radio and gamma-ray ephemerides 
as expected, further supporting our detection and AG- 
ILE timing calibration. 

We observed PSRB1509-58 in three energy bands: 
30-100 MeV, 100-500 MeV and above 500 MeV. 
We did not detect pulsed emission at a significance 
CT > 2 for £; > 500 MeV. The 7-ray light curves of 
PSRB1509-58 for different energy bands are shown 
in Fig. [H The AGILE E > 30 MeV light-curve 
shows two peaks at phases (pi — 0.39 ± 0.02 and 
(j)2 = 0.94 ± 0.03 with respect to the single radio peak, 
here put at phase 0. The phases are calculated using a 
Gaussian fit to the peaks, yielding a FWHM of 0.29(6) 
for the first peak and of 0.13(7) for the second peak. 
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FIG. 2: Fermi-LAT phase-aligned gamma-ray light curve 
of PSRB1509-58 at energies < 150 MeV, with radio peak 
at phase 0. 



where we quote in parentheses (here and throughout 
the paper) the Ict error on the last digit. The first 
peak is coincident in phase with COMPTEUs peak 
[2"3l|. In its highest energy band (10-30 MeV) COMP- 
TEL showed the indication of a second peak (even 
though the modulation had low significance, 2. la). 
This second peak is coincident in phase with AGILE'S 
second peak (Fig. [T]). AGILE thus confirms the pre- 
viously marginal detection of a second peak. 

We also analyzed Fermi-LAT data in the same in- 
terval as covered by the AGILE observations and the 
radio ephemeris: from the beginning of the mission 
until October 2009 (longer than in the published anal- 
ysis from the Fermi-LAT Collaboration, [l]). Fermi- 
LAT data were extracted through the Science Sup- 
port Center[5l| from a region of interest (ROI) of 15 
degrees, but only a 5 degrees extraction radius was 
used for the timing analysis. We selected the events 
and the correct time intervals using the Science Tools. 
We used the "diffuse" class events (highest probabil- 
ity of being gamma-ray photons) under the P6_V3 
instrument response function (IRFs), and excluded 
events with zenith angles > 105° to reject atmospheric 
gamma-rays from the Earths limb. The events were 
selected using the standard software package Science 
Tools-09-21-00 for the Fermz-LAT data analysis: in 
particular gtselect and gtmktime for the selection of 
photons and time intervals, and gthary to barycenter 
the photons. The actual photon folding is done using 
the AGILE software for pulsar observations. The re- 
sulting light curve in the soft energy band [E < 150 
MeV) are presented in Fig. [5] and it clearly shows two 
peaks. 

Fig. El shows the SED of PSRB1509-58 based on 
AGILE'S and GOMPTEUs observed fluxes. Fermi- 
LAT upper limits from [Ij] are also shown, which are 
consistent with our measurements at a 2ct confidence 
level. GOMPTEL observed this pulsar in three energy 
bands: 0.75-3 MeV, 3-10 MeV, 10-30 MeV, suggest- 
ing a spectral break between 10 and 30 MeV. AGILE 
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FIG. 3: SED of PSRB1509-58 (solid line) obtained from 
a fit of pulsed fluxes from soft to hard gamma-rays. The 
three round points represent COMPTEL observations [2^ . 
The two square points represent A GILE pulsed flux in two 
bands (30 < S < 100 MeV and IQO < E < 500 MeV). 
The thicker arrow represents AGILE upper limit above 
500 MeV. The two thin arrows represent Fermi upper lim- 
its from [l| 

pulsed flux confirms the presence of a soft spectral 
break. As shown in Fig. [31 we modeled the observed 
COMPTEL and AGILE fluxes with a power- law 
plus cutoff fit using the Minuit minimization package 
(James et al. 1975): F{E) = kx E-"" eyi^[-{E / E^Y], 
with three free parameters: the normalization fc, the 
spectral index a, the cutoff energy Ec and allowing 
/3 to assume values of 1 and 2 (indicating either an 
exponential or a superexponential cutoff). No accept- 
able values were obtained for a superexponential 
cutoff, the presence of which can be excluded at a 
3.5(7 confidence level, while for an exponential cutoff 
we found xt = 3.2 for = 2 degrees of freedom, corre- 
sponding to a null hypothesis probability of 0.05. The 
best values thus obtained for the parameters of the 
fit were: k = 1.0(2) x 10"^ cm-^, a ^ 1.87(9), 
E^ = 81(20) MeV. 

III. DISCUSSION 

The bulk of the spin-powered pulsar flux is usually 
emitted in the MeV-GeV energy band with spectral 
breaks at < 10 GeV (e.g. (2). PSRB1509-58 has the 
softest spectrum observed among gamma-ray pulsars, 
with a sub-GeV cutoff at £; w 80 MeV. In the fol- 
lowing we discuss how the new AGILE observations 
can constrain the models for emission from the puls ar 
magnetosphere (for an extended discussion see [35|)- 

When PSRB1509-58 was detected in soft gamma- 
rays but not signiflcantly at i? > 30 MeV, it was pro- 
posed that the mechanism responsible for this low- 
energy spectral break niight be photon splitting [l3|. 
The photon splitting [3| is a third-order quantum 
electro-dynamics process expected when the magnetic 
field approaches or exceeds the critical value defined 



as B„ = mlc^/{eh) = 4.413 x lO^^ q ^ost cur- 
rent theories for the generation of coherent radio emis- 
sion in pulsar magnetospheres require formation of an 
electron-positron pair plasma developing via electro- 
magnetic cascades. In very high magnetic fields the 
formation of pair cascades can be altered by the pro- 
cess of photon splitting: 7 — > 77, which will oper- 
ate as an attenuation mechanism in the high-field re- 
gions near pulsar polar caps. Since it has no energy 
threshold, photon splitting can attenuate photons be- 
low the threshold for pair production, thus determin- 
ing a spectral cutoff at lower energies. 

In the case of PSRB1509~58 a polar cap model 
with photon splitting would be able to explain the 
soft gamma-ray emission and the low energy spectral 
cutoff, now quantifled by AGILE observations. Based 
on the observed cutoff, which is related to the photons' 
saturation escape energy, we can derive constraints on 
the magnetic field strength at emission, in the frame- 
work of photon splitting: 

e^:*~O.O77(B'sin0fcB,o)"'/' (1) 

where e^"* is the photon saturation escape energy, 
B' = B/Bcr and 9kB,o is the angle between the 
photon momentum and the magnetic field vectors at 
the surface and is here assumed to be very small: 
dkB,o < 0.57° (see [13 )• Using the observed energy 
cutoff {elf^ ~ = 80 MeV) we find that B' > 0.3, 
which implies an emission altitude < 1.3Rns, which 
is the height where possibly also pair production could 
ensue. This altitude of emission agrees with the polar 
cap models (see e.g. [l0|)- A smaller energy cutoff, 
as in (17| . would have implied even lower emission al- 
titude and a sharper break, possibly caused by the 
total absence of pair production. It is apparent that 
small differences in the emission position will cause 
strong differences in spectral shape. This is possi- 
bly the reason for the different emission properties of 
the two peaks as observed in the total {AGILE plus 
GOMPTEL) gamma-ray energy band and confirmed 
by our reanalysis of Fermz-LAT data where this peaks 
clearly appears in the soft energy band. Also, a trend 
can be observed, from lower to higher energies (see 
the X-ray light-curve for the trend in the first peak, 
as in Fig. 3 of [23]), of the peaks slightly drifting away 
from the radio peak. This we assume to be another 
signature of the fact that small variations in emis- 
sion height can be responsible for sensible changes in 
the light curves in such a high magnetic field. The 
scenario proposed by [l3l is strengthened by its pre- 
diction that PSR B0656-I-14 should have a cutoff with 
an intermediate value between PSRB1509-58 and the 
other gamma-ray pulsars. Additionally, PSRB1509- 
58 [1, i3 and PSR B0656-hl4 [ll|, show evidence 
of an aligned geometry, which could imply polar cap 
emission. 

The polar cap model as an emission mechanism is 
debated. From the theoretical point of view, the an- 
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gular momentum is not conserved in polar cap emis- 
sion (see [4i|| for a revision). Also, a preferential 
explanation of the observed gamma-ray light curves 
with high altitude cascades comes from the recent re- 
sults by Fermi-LKT In the case of PSRB1509- 
58, the derived gamma-ray luminosity from the flux 
at i? > 1 MeV, considering a 1 sr beam sweep is 
L-y — 4.2l;Q2d§ 2 ^ 10^^ erg/s, where ^5.2 indicates 
the distance in units of 5.2 kpc. While tradition- 
ally the beaming fraction (/o) was considered to be 
the equivalent of a 1 sr sweep, nowadays (see e.g. 
(46j ) the tendency is to consider a larger beaming 
fraction (/n ss 1), close to a An sr beam. Using 
/o = 1 in our calculations, we would have obtained 
= 5.8^0 8^5 2 X 10^^ erg s^^ Thus the maximum 
conversion efficiency of the rotational energy loss {E w 
1.8 X 10'^^ erg s~^) into gamma-ray luminosity is 0.3. 
Our result is not easily comparable with the typical 
gamma-ray luminosities above 100 MeV, because for 
PSRB1509-58 this energy band is beyond the spec- 
tral break. Using AGILE data alone we obtained a 
luminosity above 30 MeV = 5.2(6)di 2 x 10^^ erg/s, 
again for a 1 sr beam. If the gamma-ray luminosity 
cannot account for a large fraction of the rotational 
energy loss, then the angular momentum conservation 
objection becomes less cogent for this pulsar, exactly 
as it happens for the radio emission. 

Alternatively, if such an efficiency as that of 
PSRB1509-58 were incompatible with this conserva- 
tion law, an interpretation of PSRB1509-58 emission 
can be sought in the frame of the three dimensional 
outer gap model [13] ■ According to their model, hard 
X-rays and low energy gamma-rays are both produced 
by synchrotron self-Compton radiation of secondary 
e+e" pairs of the outer gap. Therefore, as observed, 
the phase offset of hard X-rays and low energy gamma- 
rays with respect to the radio pulse is the same, with 
the possibility of a small lag due to the thickness of 
the emission region. According to their estimates a 
magnetic inclination angle a ~ 60° and a viewing an- 
gle C « 75° are required to reproduce the observed 
light curve. Finally, using the simulations of (46l | , who 
produced a map of pulse profiles for different com- 
binations of angles a and C, in the different emission 
models, the observed light curve from AGILE is best 
reproduced if a ~ 35° and ^ ~ 90°, in the framework 
of the two pole caustic model [isj . 



The values of a and C, required by the [50[ model are 
not in good agreement with the corresponding values 
obtained with radio measurements. In fact, observe 
that a must be < 60° at the 3(T level. The prediction 
obtained by the simulations of (46ll better agrees with 
the radio polarization observations. In fact, in the 
framework of the rotating vector model (RVM, see 



e.g. [25| and references therein), [9j also propose that, 
if the restriction is imposed that C, > 70° [2^, then 
a > 30° at the 3a level. For these values, however, 
the Melatos model for the spin down of an oblique 
rotator predicts a braking index n > 2.86, slightly 
inconsistent with the observed value {n = 2.839(3), 
see [ii]). Also in the case of PSR B0656-I-14, [4| 
conclude that the large values of a and C are some- 
what at odds with the constraints from the model- 
ing of the radio data and the thermal X-rays which 
seem to imply a more aligned geometry. Improved 
radio polarization measurements would help placing 
better constraints on the pulsar geometry and there- 
fore on the possibility of a gap in the extended or 
outer magnetosphere, but the quality of the polariza- 
tion measurements from is already excellent, the 
problem being that PSRB1509-58, like most pulsars, 
only shows emission over a limited pule phase range 
and therefore the RVM models are highly degenerate. 

At present the geometry privileged by the state of 
the art measurements is best compatible with polar 
cap models. Higher statistics in the number of ob- 
served gamma-ray pulsars could help characterize a 
class of "outliers" having gamma-ray emission from 
the polar caps, which potentially constitute a privi- 
leged target for AGILE. 
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